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Topological Shocks in Burgers Turbulence 
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|V^| 2 = !/V 2 ^ + F(i,t). 
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The dynamics of the multi-dimensional randomly forced Burgers equation is studied in the limit of 
vanishing viscosity. It is shown both theoretically and numerically that the shocks have a universal 
global structure which is determined by the topology of the configuration space. This structure is 
shown to be particularly rigid for the case of periodic boundary conditions. 
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The d— dimensional randomly forced Burgers equation 

d t u+(u- V)u = vX7 2 u- VF(x,t), (1) 

appears in a number of physical problems, ranging from 
the dynamics of interfaces and cosmology to hydrody- 
namics (see, e.g., Ref. jl| for a review). In the context 
of fluid dynamics, the Burgers equation, frequently re- 
ferred to as "Burgers turbulence" , is a simple model for 
analyzing the signature of singularities, mostly shock dis- 
continuities, in the statistics of the velocity field (see, e.g., 
Refs. [§, |, || §, §]). The (statistical) steady-state theory 
for Burgers turbulence in the limit of vanishing viscos- 
ity (y — * 0) was developed for d = 1 in the spatially 
periodic case . The analysis of the Lagrangian dynam- 
ics led to the distinguishing of a particular trajectory, 
the global minimizer, corresponding to the unique fluid 
particle that is never absorbed by a shock. The counter- 
part to the global minimizer is a unique main shock, with 
which all other shocks arc merging after a finite time and 
hence, in which all the matter gets concentrated. Here 
the goal is to show that these objects, extended to multi- 
dimensional situations, determine the global structure of 
the stationary solution. This structure is strongly con- 
nected with the topology of the configuration manifold 
defined by the boundary conditions. We show that in 
any dimension, a unique global minimizer exists, so that 
the shocks have either a local or a global topological na- 
ture. The global shocks, unavoidably present in Burgers 
dynamics, are called the topological shocks; they have a 
nontrivial structure for d > 1. For simplicity, we mostly 
consider space-periodic forcing potentials for which the 
configuration manifold is the ci-dimensional torus T d (1- 
periodic boundary conditions), but most of our work can 
be extended to other types of boundary conditions and 
configuration spaces. 

If the initial data is of gradient type, the velocity 
field preserves this property at any later time, so that 
u(x,t) = —Vtp(x,t), where the velocity potential ip 



(2) 



When the external potential F is delta-correlated in both 
space and time, this equation is known as the Kardar- 
Parisi-Zhang model for interface dynamics Jjj . We focus 
here on smooth-in-space forcing potentials with a corre- 
lation function given by 

(F(x 1 ,t 1 )F(x 2 ,t 2 )) = G{ Xl - x 2 ) 5{h - i 2 ), (3) 

where G is a smooth large-scale function. This type of 
large-scale forcing was chosen analogous with that usu- 
ally assumed in work on forced Navier-Stokes turbulence. 
Note that the results discussed in this Letter can be ex- 
tended to other types of random forcing (e.g. with a finite 
correlation time). Because of space periodicity, the av- 
erage velocity b = J Jd u(x, t)dx is an integral of motion. 
Its value does not affect the structure of the topological 
shocks and, for simplicity, we choose 6 = 0. 

The initial-value problem associated to the Hamilton- 
Jacobi equation (|^), in the inviscid limit v — > 0, has a 
variational solution || . Denoting tpo the potential at the 
initial time to, the velocity potential at times t > to is 
given by 



ip{x,t) 
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where the infimum is taken over all differentiable curves 
7 : [to,t] — > T d such that *y(t) = x and A is the La- 
grangian action 



A(j,t) 



la 



i| 7 (r)| 2 -F( 7 (r),r 



dr. 



(5) 



A minimizing trajectory is called a minimizer on the in- 
terval [to,t] and is associated to a fluid particle reaching 
x at time t. In the limit to — > -co, a stationary regime 
is obtained, independent of i/V The solution is then de- 
termined by one-sided minimizers, i.e. action-minimizing 
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trajectories from — oo to t. It is easily seen from Eq. 
(0) that all the minimizers are solutions of the Euler- 
Lagrange equations 



7(r) = v(t), 

v(t) = -VF( 7 (t),t) 



(6) 
(7) 



A global minimizer (or a two-sided minimizer) is defined 
as a curve which minimizes the action for any time in- 
terval [ii,i2] and thus corresponds to the trajectory of a 
fluid particle that is never absorbed by shocks. 

We now state three main results whose rigorous proof 
is given in Ref. || and that are essential for the intro- 
duction of topological shocks. First, there exists a unique 
solution of the Hamilton-Jacobi equation (^) in the limit 
v — > which is extendible to all times. This solution is 
continuous and almost everywhere differentiable in space. 
It generates uniquely a stationary distribution for the 
random Hamilton-Jacobi equation and its gradient de- 
fines a unique statistically stationary solution of the in- 
viscid Burgers equation. Second, for a given time and 
for every space location where the potential is differen- 
tiable, there exists a unique one-sided minimizer. The 
locations where the one-sided minimizers are not unique 
correspond to shock positions. Finally, there exists a 
unique global minimizer. This third statement is cru- 
cial for the construction of topological shocks because it 
implies that for large negative times r — > — oo, all the 
one-sided minimizers are asymptotic to the global mini- 
mizer po[ . 

To introduce the notion of topological shock, we "un- 
wrap" , at a given time t, the configuration space T 2 to 



l(af ). Now, for a given 



the entire space R (see Fig. 
realization of the forcing, we obtain instead of a single 
global minimizer an infinite number of them, each being 
the image of others by integer shifts. They form a lattice 
parameterized by vectors k with integer components and 
are denoted 7^ g '. The backward- in-time convergence to 
the global minimizer on T d implies that every one-sided 
minimizer emanating from some location in R d is asymp- 
totic to a particular global minimizer "/^P on the lattice. 
Hence, every location x which has a unique one-sided 
minimizer is associated to an integer vector k(x), defin- 
ing a tiling of the space at time t. The tiles Ok are 
the sets of points whose associated one-sided minimizer 
is asymptotic to the fc-th global minimizer. The bound- 
aries of the Ofc's correspond to the positions of particular 
shocks that are called the topological shocks. They are 
the locations for which at least two one-sided minimizers 
approach different global minimizers on the lattice. In- 
deed, a point where two tiles Ok t and Ok 2 meet, has at 
least two one-sided minimizers, one of which is asymp- 
totic to 7^ and another to 7^ • Of course, there are 
also points on the boundaries where three or more tiles 
meet and thus where more than two one-sided minimiz- 
ers are asymptotic to different global minimizers. For 



d = 2 such locations are isolated points corresponding 
to the intersections of three or more topological shock 
lines, while for d = 3, they form edges and vertices where 
shock surfaces meet. Note that, generically, there exist 
other points inside Ok with several minimizers. They 
correspond to shocks of "local" nature because at these 
locations, all the one-sided minimizers are asymptotic to 
the same global minimizer 7^ and hence, to each other. 
In terms of mass dynamics, the topological shocks play 
a role dual to that of the global minimizer. Indeed, all 
the fluid particles are converging backward-in-time to the 
global minimizer and are absorbed forward-in-time by 
the topological shocks. Assuming that the Burgers equa- 
tion (0) is accompanied by a continuity equation for the 
mass density, this implies that all the mass concentrate 
at large times in the topological shocks. 
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FIG. 1: (a) Sketch in space-time of the unwrapped picture 
for d = 2; on the horizontal plane which corresponds to a 
lattice of periodic box replicas at time t, the bold lines denote 
the topological shocks and triple points are represented by 
dots. The different "horns" illustrate the backward-in-time 
convergence to four different global minimizers represented 
as dashed lines and the filled areas represent different tiles 
Ok- (b) Position of the topological shock on the torus; the 
two triple points are represented as dots. 

We now describe the global structure of the topologi- 
cal shocks. Parameter counting suggests that there are 
generically (d— l)-dimensional surfaces of points with two 
one-sided minimizers which contain (d — 2)-dimensional 
sub-manifolds with three one-sided minimizers and so on. 
This ends up with single points (zero-dimension) from 
which emanate (d + 1) one-sided minimizers. As one ex- 
pects to see only generic behavior in a random situation, 
the probability to have points associated to more than 
(d+l) one-sided minimizers is zero. It follows that there 
are no points where (d + 2) tiles Ok meet, an important 
restriction on the structure of the tiling. Thus for d = 2, 
the tiling is constituted of curvilinear hexagons. Indeed, 
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suppose each tile Ok is a curvilinear polygon with s ver- 
tices corresponding to triple points. For a large piece 
of the tiling which consists of N tiles, the total num- 
ber of vertices is n v ~ siV/3 and the total number of 
edges is n e ~ sN/2. The Euler formula implies that 
1 = n v — n e + N ~ (6 — s)N/6; so we have s = 6, corre- 
sponding to an hexagonal tiling. As shown in Fig. l(b% 
this structure corresponds on the periodicity torus T 2 , 
to two triple points connected by three shock lines which 
are the curvilinear edges of the hexagon Oq. The connec- 
tion between the steady state velocity potential and the 
topological shocks is shown on Fig. ^ which was obtained 
numerically. 




FIG. 2: Snapshot of the velocity potential ip(x,y,i) for d — 2 
in the statistical steady state, obtained numerically with 256 2 
grid points. Shock lines, corresponding to locations where ip is 
not differentiable, are represented as black lines on the bottom 
of the picture; the four gray areas are different tiles separated 
by the topological shocks; the other lines are local shocks. 

Although the topological shocks always form hexago- 
nal patterns when d — 2, the corresponding tilings can be 
of different types; in the course of time, the merger of two 
triple points is the generic mechanism for changing the 
type of the tiling. This so-called flipping bifurcation H] 
has the property of redistributing matter among nodes, 
so that the mass does not concentrate in a particular 
triple point. In higher dimensions, the structure of topo- 
logical shocks can be more complex; for instance, it is not 
possible to determine uniquely the shape of polyhedra 
forming the tiling. Nevertheless, the minimal polyhedra 
defining such tilings for d = 3 can be shown to have 24 
vertices @ 

All the above results concerning the global structure 
of solutions require a statistical steady state, achieved 
asymptotically at large times. The convergence to this 
regime is actually exponential so that, generally, the 
global picture of the flow is reached after just a few 
turnover times. The nature of the convergence is related 
to the local properties of the global minimizer and more 
particularly to its hyperbolicity. For d = 1, the global 
minimizer has been shown to be a hyperbolic trajectory 



of the dynamical system defined by the Euler-Lagrange 
equations (§-0) 0. In multi-dimensional situations the 
hyperbolicity of the global minimizer is an open problem. 
Since the Lagrangian flow defined by (^-0) is Hamilto- 
nian, one can define d pairs of non-random Lyapunov ex- 
ponents with opposite signs. Hyperbolicity means that 
none of these exponents vanishes. This question can be 
addressed in terms of the backward-in-time convergence 
of the one-sided minimizers to the global one or, in terms 
of forward-in-time dynamics, by looking at how fast La- 
grangian fluid particles are absorbed into shocks. For 
this, we consider the set f2(T) of locations x such that 
the fluid particle at x at time t — survives, i.e. is not 
absorbed by any shock, until the time t = T. The long- 
time shrinking of f2 as a function of time is asymptot- 
ically governed by the Lyapunov exponents. To ensure 
the absence of vanishing Lyapunov exponents, it is suf- 
ficient to show that the diameter of Q(T) decays expo- 
nentially as T — > oo. Below we demonstrate numerically 
that this is indeed the case for d = 2. For this we assume 
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FIG. 3: Time evolution of the diameter of the Lagrangian set 
f2(T) (points corresponding to the regular region) for three 
different types of forcing spectra normalized to give the same 
Lyapunov exponents; average over 100 realizations and with 
256 2 grid points. 



that the forcing is a sum of independent random impulses 
concentrated at discrete times a case to which the 
present theory remains applicable. Between "kicks" the 
velocity field decays according to the unforced Burgers 
equation. This allows us to use the fast Legendre trans- 
form method |l5[ , based on discrete approximations of 
the one-sided minimizers over a grid, which gives directly 
the solution in the inviscid limit and is particularly useful 
due to its strong connection with the Lagrangian picture 
of the flow. We can then track numerically the set f2(T) 
of regular Lagrangian locations. As shown in Fig. for 
three different types of forcings, the diameter of this set 
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decays exponentially in time, providing strong evidence 
for the hyperbolicity of the global minimizer when d = 2. 

Since all the one-sided minimizers converge backward- 
in-time to the global minimizer, hyperbolicity implies 
that, in the statistical steady state, the graph of the 
solution in the phase space (as, u) is made of pieces of 
the smooth unstable manifold associated to the global 
minimizer with discontinuities along the shocks lines or 
surfaces. In other words, shocks appear as jumps be- 
tween two different folds of the unstable manifold. The 
smoothness of the unstable manifold is key; for instance, 
it implies that when d = 2, the topological shock lines are 
smooth curves. The above geometrical construction of 
the solution has much in common with that appearing in 
the unforced problem. Indeed, when F = 0, the solution 
can be obtained by considering in the (as, it) space, the 
Lagrangian manifold defined by the position and the ve- 
locity of the fluid particles at a given time. This analogy 
gives good ground to conjecturing that several universal 
properties associated to the unforced problem still hold 
in the forced indeed happens in one dimension 

j|, 0). Hyperbolicity implies that there exists a strong 
parallel between the forced and the unforced situations. 
Despite the fact that the mass dynamics is completely 
different in the two cases (mass is absorbed by shocks lin- 
early in time in the decaying case and exponentially fast 
in the forced case) , many features of the velocity field are 
universal. For instance, it was shown in Ref. |l(| that, 
for the unforced case and d > 1, large but finite mass 
densities are localized near time-persistent boundaries of 
shocks ( "kurtoparabolic" singularities) contributing, in 
any dimension, power-law tails with the exponent — 7/2 
in the probability density function (PDF) of both veloc- 
ity gradients and mass densities. When a force is applied, 
the geometry of the solutions is very similar to that ap- 
pearing in unforced situations. This leads again to a 
universal —7/2 power-law behavior of the PDF of veloc- 
ity gradients and mass density, irrespective of d. The 
issue of similarities between forced and unforced Navier- 
Stokes turbulence and the search for universal statistical 
properties of the velocity field is, of course, still an open 
problem. 

Note, finally, an important physical problem which will 
be addressed in a future work, namely the understand- 
ing of the behavior of minimizers and the effects related 
to global shocks in the case of spatially extended non- 
periodic systems. This amounts to investigating inter- 
mediate time asymptotics when the size of the system is 
much larger than the forcing scale. Preliminary numer- 
ical results indicate the appearance of a shock structure 
resembling the structure of topological shocks. 
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